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Abstract 

A brief overview is given over recent results on the spectral properties of spherically symmetric MHD 
a 2 — dynamos. In particular, the spectra of sphere-confined fluid or plasma configurations with physically 
realistic boundary conditions (BCs) (surrounding vacuum) and with idealized BCs (super-conducting sur- 
rounding) are discussed. The subjects comprise third-order branch points of the spectrum, self-adjointness 
of the dynamo operator in a Krein space as well as the resonant unfolding of diabolical points. It is sketched 
how certain classes of dynamos with a strongly localized a— profile embedded in a conducting surrounding 
can be mode decoupled by a diagonalization of the dynamo operator matrix. A mapping of the dynamo 
eigenvalue problem to that of a quantum mechanical Hamiltonian with energy dependent potential is used to 
obtain qualitative information about the spectral behavior. Links to supersymmetric Quantum Mechanics 
and to the Dirac equation are indicated. 

Preliminaries 

The magnetic fields of planets, stars and galaxies are maintained by dynamo effects in conducting fluids or 
plasmas [TJ [21 [3] . These dynamo effects are caused by a topologically nontrivial interplay of fluid (plasma) 
motions and a balanced self- amplification of the magnetic fields — and can be described within the framework 
of magnetohydrodynamics (MHD) [TJ [5] . 

For physically realistic dynamos the coupled system of Maxwell and Navier-Stokes equations has, in general, 
to be solved numerically. For a qualitative understanding of the occurring effects semi-analytically solvable toy 
models play an important role. One of the simplest dynamo models is the so called a 2 — dynamo with spherically 
symmetric a— profilcQ a(r) (see, e.g. [2]). For such a dynamo the magnetic fields can be decomposed in poloidal 
and toroidal components, expanded over spherical harmonics 011] and unitarily re-scaled [5]. As result one 
arrives at a set of /—mode decoupled 2 x 2— matrix differential eigenvalue problems [3J[U[5] 

«a = ( _q {1] ) , Q[a] ■= -d r a{r)d r + a(r)^±ll (i) 

with boundary conditions (BCs) which have to be imposed in dependence of the concrete physical setup and 
which will be discussed below. The a— profile describes the net effect of small scale helical turbulence on the 
magnetic field [2]. It can be assumed real-valued a(r) £ K and sufficiently smooth. We note that the reality 
of the differential expression ([1]), independently from the concrete BCs, implies an operator spectrum which is 
symmetric with regard to the real axis, i.e. which consists of purely real eigenvalues and of complex conjugate 
eigenvalue pairs. 
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1 The a— profile a(r) plays the role of an effective potential for the a 2 — dynamo. 
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In [1] it was shown that the differential expression ([TJ) of this operator has the fundamental (canonical) 
symmetry [BJ[7] 

'-(1 '»)■ P> 



In case of BCs compatible with this fundamental symmetry the operator turns out self-adjoint in a Krein spacc0 
(ICj, [., [HIS] and in this way it behaves similar like Hamiltonians of VT— symmetric Quantum Mechanics 

(ptsqm) ma nn Ha na na . 
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Figure 1: Real and imaginary components of the a 2 — dynamo spectrum as functions of the scale factor C of an 
a— profile a(r) = C x (1 — 26.09 x r 2 + 53.64 x r 3 — 28.22 x r 4 ) in the case of angular mode number I = 1 and 
physically realistic boundary conditions (J3J). The concrete coefficients in the quartic polynomial a(r) have their 
origin in numerical simulations of field reversal dynamics (see Ref. [201 121)). Only the imaginary components 
with > are shown. The symmetrically located complex conjugate (3A < 0)— components are omitted for 
sake of brevity. 



Subsequently, we first present a sketchy overview of some recent results on the spectral behavior of a 2 — dynamos 
obtained in [SJ HH HH1 EZ1 HS1 US] which we extend by a discussion of the transition from a 2 — dynamo configu- 
rations confined in a box to dynamos living in an unconfined conducting surrounding. 
Physically realistic BCs and spectral triple points 

For roughly spherically symmetric dynamical systems like the Earth the conducting fluid is necessarily confined 
within the core of the Earth so that the a— effect resulting from the fluid motion has to be confined to this 
core. Setting the surface of the outer core at a radius r = 1 one can assume a(r > 1) = and a behavior of the 
magnetic field at r > 1 like in vacuum. A multi-pole-like decay of the magnetic field at r — > oo leads then to 

2 For comprehensive discussions of operators in Krein spaces see, e.g., [6]0[8]. 
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mixed effective BCs at r = 1 (see, e.g. [2]) and an corresponding operator domain of the type 
2?(2t Q ) - {u e H = L 2 (0, 1) © L 2 (0, 1)| u(r \ 0) = 0, <Bu| r ^i = o} , 
35 := 



u := 



"l 

"2 



0r+ l - 

1 



(3) 



From the domain 2?(2lJ,) of the adjoint operator 

2?(»+ ) = {ue?i = L a (0, 1) © L a (0, 1)1 u(r \ 0) = 0, «8u| r/ >i = o} , 
03 := 



u := 



U-2 



d r + i -a(r)9 r 
1 



(4) 



one reads off that 2?(2ljj ^ 2?(2l Q ) and, hence, the dynamo operator % a itself is not self-adjoint even in a Krein 
space. 
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J -dynamo with a(r) = C [-(21.465 + 2.467C) + (426.412 + 167.928<> 2 - (806.729 + 436.289C)r 3 
272.991£)r 4 ] and a spectral triple point at (£ = 0.45, C — 0.86). Highlighted (fat) lines correspond 



Figure 2: n 
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to purely real branches of the spectrum. The cusp in the imaginary component (lower right graphics) indicates 
the closely located triple point 



In case of constant a— profiles and arbitrary I £ N, the spectrum is implicitly given by a characteristic 
equation built from spherical Bessel functions [2]. In all other cases numerical studies are required. A typical 
spectral branch graph is depicted in Fig. [TJ Obviously, for the specific a— profile it contains a large number of 
spectral phase transitions from real spectral branches to complex ones and back. There are strong indications 
that phase transition points (second order branch points/exceptional points) of the spectrum close to the A = 
line play an important role in polarity reversals of the magnetic field (see [201 1211 12H H3] for numerical studies 
and [Mj for recent experiments). 

Apart from the second-order branch points visible in Fig. [T] there may occur third- and higher-order branch 
points. They are located on hyper-surfaces of higher co-dimension in parameter space and they therefore require 
a tuning of more parameters to pin them dowr0. Corresponding results have been obtained in [16] and are 
illustrated in Fig. [2] The triple points result from coalescing second-order branch points, correspond to 3 x 3 
Jordan blocks in the spectral decomposition of the operator and are accompanied by a merging or disconnecting 
of two complex spectral sectors over the parameter space. An implicit indication of a closely located triple point 
is the presence of cusps in the imaginary components as they are visible in Figs. [TJ[2] 



3 An explicit hyper-surface parametrization of second-order branch point configurations embedded in a VT— symmetric 3 X 
3— matrix model with corresponding 2 X 2— Jordan-block preserving modes can be found e.g. in the recent work 1251 . 
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Idealized BCs and Krein-space related perturbation theory 

In order to gain some deeper insight into possible dynamo-related processes semi-analytical toy model consid- 
erations play a crucial role. A certain simplification of the eigenvalue problem has been achieved in [171 118] by 
considering a reduced and idealized (auxiliary) problerrQ with Dirichlet BCs imposed at r = 1, i.e. by setting 
u(r = 1) = 0. In this case it holds 2?(2l Q ) = 2?(2lJ,) and the operator 2l Q is self-adjoint in a Krein space 
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Figure 3: The spectral mesh of the operator matrix 2l Q for 1 = (a); its resonant deformation due to harmonic 
perturbations of a constant a-profile (b), (c); the formation of overcritical oscillatory dynamo regimes for e 
increasing from e = (QA = 0) to e = 35 (5RA > 0, QA ^ for some branches) (d); and the resonant unfolding 
of DPs in the complex plane (e), (f). 

(Xj,[.,.]j) [5]. For constant a— profiles a(r) = ao = const the eigenvalue problem (2l ao — A)u = becomes 
exactly solvable in terms of orthonomalized Riccati-Bessel functions 

u n {r) = N n r 1/2 J l+ i (Jp^r), N n := (u m ,u n ) = S mn , \\u„\\ = 1 (5) 

with p n > the squares of Bessel function roots J /+ i (y 7 /^) = 0. The solutions of the eigenvalue problem have 
the form 

"""(iV^)"™ 682 ®^ ' 1 )' (6) 

are Krein space orthonormalized 

k,"n] = ± 2 ^mn, [u^,U^]=0, U± € K± C K., =: U% £ = ± (7) 

and correspond to eigenvalue branches A^ = — p n + sao^/p^ which scale linearly with ao- In the (ao, 5RA)— plane 
the branches A+ and A~ of states u+, u~ of positive and negative Krein space type form a spectral mesh 
(see Fig. |3|). The intersection points (nodes of the mesh) are semisimple double eigenvalues, i.e. eigenvalues 
of geometrical and algebraical multiplicity two — so called diabolical points (DPs) [26 . Two given branches 
A^(ao) and A^(ao) intersect at the single point A = Aq := e8-Jp n p m , = :— £y/p^+Sy/p^ and one obtains 
that branches from states of opposite Krein space type e = — 6 intersect for Ag < 0, whereas states of the same 
type (e = 6) intersect at Aq > 0. Under small inhomogeneous perturbations a(r) = ajJ + Aa(r) = cto+e<fi(r) the 

4 From a physical point of view such a 2 — dynamos can be regarded as embedded in a superconducting surrounding. 
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diabolical points split Aq h- > Ag+eAi + . . . into two real or complex points (see also [S7J for similar considerations) 
with leading contribution Ai resulting from the quadratic equation 

A 2 Ai (e [<B< ' <] I 5 [<BU "" U " ] 1 I c5 [<Bu "' u " ] [<Bu "' ~ [<Bu "' u " ] ' - Q (8) 



where 



= / v 

Jo 



s5y/p n p m H 5 — 



dr. (9) 



The unfolding of the DPs follows the typical Krein space rule. When they result from branches of the same 
type (e — 5,Xq > 0) then the corresponding DPs unfold purely real-valued, whereas DPs from branches of 
opposite type (e = — <5, Aq < 0) may unfold into complex conjugate eigenvalue pairs. This behavior is clearly 
visible in Fig. [3] b,c. Direct inspection reveals that the spectral meshes of unperturbed operators 2t Qo for 
I = and < I <C oo show strong qualitative similarities so that results obtained for the quasi-exactly solvable 
(I = 0)— model will qualitatively hold for models with < I -C oo too. Via Fourier expansion of a(r) a very 
pronounced resonance has been found along parabolas in the (ao, 5RA)— plane indicated by white and colored dots 
in Fig. — leaving regions away from these parabolas almost unaffected. An especially pronounced resonance 
is induced by cosine perturbations which in linear approximation affect only the single parabola j = 2k, Fig. 
[5jD,e. Sine perturbations act strongest on parabolas |j| = 2k ± 1 with decreasing effect on \j\ — 2k ± m for 
increasing m (see Fig. [3J;,f). Physically important is the fact that higher mode numbers k (shorter wave 
lengths of the Aa(r) perturbations) affect more negative 5RA. Due to a magnetic field behavior cx e xt this is 
the mathematical formulation of the physically plausible fact that small-scale perturbations decay faster than 
large-scale perturbations. Numerical indications for the importance of this behavior in the subtle interplay of 
polarity reversals and so called excursions (" aborted" reversals) of the magnetic field have been recently given 
in [23]. 

Diagonalizable a 2 dynamo operators, SUSYQM and the Dirac equation 

Another approach to obtain quasi-exact solution classes of the eigenvalue problem (2l a — A)u = consists in a 
A— dependent diagonalization of the operator matrix ([T]). The basic feature of this technique, as demonstrated in 
19J, is a two-step procedure consisting of a gauge transformation which diagonalizes the kinetic term and a sub- 
sequent global (coordinate-independent) diagonalization of the potential term. Such an operator diagonalization 
is possible for a— profiles satisfying the constraint 

a"(r) + ^a 3 (r) - a 2 a(r) = (10) 

with a = const G R a free parameter. Solutions a(r) of this autonomous differential equation (DE) can be 
expressed in terms of elliptic integrals. In order to maximally explore similarities to known QM type models^ 
a strongly localized a— profile has been assumed which smoothly vanishes toward r — > oo. Physically, such a 
setup can be imagined as a strongly localized dynamo-maintaining turbulent fluid/plasma motion embedded 
in an unbounded conducting surrounding (plasma) with fixed homogeneous conductivity. The only a— profile 
with a(r — + oo) — * satisfying (|10[) has the form of a Korteweg-de Vries(KdV)-type one-soliton potential 

a M = ^T7^ vi ■ ( n ) 

cosh|a(r — ro)J 

This amazing finding indicates on deep structural links to KdV and supersymmctric quantum mechanics 
(SUSYQM) and opens up a completely new exploration approach to a 2 — dynamos^. In [19] we restricted 
the consideration to the most elementary solution properties of such models. The decoupled equation set after 
a parameter and coordinate rescaling has been found in terms of two quadratic pencils 

[-dl + l -^l- 1 -a 2 + 1 - T ea-e 2 }F ± =0, a = (12) 

or 2 2 cosh(x — Xo) 

1 /2 

in the new variable x :— ar and with new auxiliary spectral parameter e = — A) . The equivalence 
transformation from (2t Q — A)u = to (fT2"]l is regular for e 7^ and becomes singular at e = where (TT2"]) has to 
be replaced by a Jordan type equation system 



di - V -Vi \ ( Ei 
d 2 x - V 



= (13) 



5 For early comments on structural links between MHD dynamo models and QM-related eigenvalue problems see e.g. 1281 . 
6 The question of whether this new class of quasi-exactly solvable a 2 — dynamo models might be structurally related (via dynamical 
embedding) to the recently studied VT— symmetrically extended KdV solitons [29| 130] remains to be clarified. 
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with potentials Vq — 1(1 + l)x 2 
eigenvalue problems (JTSJ) read 



i(« 2 - 1), Vi 



In terms of the original spectral parameter A the 



1(1 + 1) 



-2 a T 



1/2 



F ± = -XF± 



(14) 



and can be related to the spectral problem of a QM Hamiltonian with energy E = — A and energy-dependent 

potential component + (\ — A) 1 ^ 2 a(x) = =p (E + |) a(a;). For physical reasons asymptotically vanishing 
field configurations with F±(x — > 00) — > 0, 30,1(0; — > 00) — > are of interest. These Dirichlet BCs at infinity 
imply the self-adjointness of the operator 2l Q in a Krein space JCj — with (| 1 2[) . (p~3|) as special representation 
of the eigenvalue problem (2l Q — A)u = 0. From the structure of (jT2|) . (fl4|) follows that the only free parameter 
apart from the angular mode number I £ N is the maximum position xq of the a— profile a(x) (the minimum 
position of the potential component —a 2 (x)/2) so that solution branches will be functions A(xq). 





Figure 4: Spectra e(xo) (a) and A(xo) (b) in case of angular mode numbers / = 0,1,2,3. For numerical reasons the 
Dirichlet BC has been imposed at the large distance x — 100. 

With the help of SUSY techniques it has been shown in 19J that (fl"4"|) has a single bound state (BS) type 
solution which via E = — X < corresponds to an overcritical dynamo mode A > 0. It has been found that the 
BS solutions of (|14[) behave differently for xq < xj and xq > xj, where for xq = xj the description in terms of 
(fT2| breaks down and has to be replaced by the singular Jordan type representation (p~3|) . By a SUSY inspired 
factorization ansatz 

„, 1(1 + 1) 1 9 1 + 
~ d l + - 2 a +2 =L L ' (15) 

L = — d x + w, L^ = d x +w, w — u'/u, (16) 

an equivalence relation between (fT!?|) and a system of two Dirac equations 

H±V± = e*± , H± = jd x + V± (17) 

— (fe)-(^). -(-1 J)" ») (18) 

has been established for models with x$ < xj. General results on Dirac equations allowed then for the conclusion 
that in case of xq < xj the bound state related spectrum has to be real. A perturbation theory with the distance 
5 = xq — xj from the Jordan configuration as small parameter supplemented by a bootstrap analysis showed 
that the Dirichlet BCs F±(x — > 00) — > render only the solution F±(x) non-trivial and with real eigenvalue, 
whereas F-(x) has to vanish identically F-(x) = 0. The single spectral branch in terms of A(xo) and e(xo) is 
depicted in Fig. [4] for angular mode numbers I = 0, 1, 2, 3. 

Assuming the dynamo model with strongly localized a— profile pip . (|12p confined in a large box, i.e. with 
Dirichlet BCs imposed at large x = X 3> 0, one can study the dynamo spectrum in the infinite box limit. 
Figures [5] a, b show the corresponding behavior. Due to its localization the BS-related overcritical dynamo mode 
is almost insensitive to the X — ► 00 limit. This is in contrast to the under-critical (decaying) modes which 
behave as expected for a sign inverted box spectrum of QM. For fixed mode number n > 2 and X — > 00 the 
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Figure 5: Cutoff (X-)dependence of the spectral branches with radial mode numbers n = 1,2 and angular mode 
numbers I = (a) and I — 2 (b) for cutoffs (box-lengths) X — 10, 20, 40. Clearly visible are the X— independence of 
the overcritical BS type modes in — 1) and the tendency A cx — 1/X 2 for the undercritical (box- type) mode. The modes 
with n > 3 show the same qualitative A oc —1/X 2 behavior like the n = 2 mode and are note depicted here. 



energies E n decrease like E n cx 1/X 2 \ and the corresponding part of the spectrum becomes quasi-continuous 
and related to the continuous (essential) spectrum of QM scattering states of a particle moving in the energy 

dependent potential '^"^ — ^a 2 (x) =p [E + ^) 1 ^ 2 a(x). For the associated dynamo eigenvalues this implies 
A„ oc —1/X 2 / — as it is clearly visible in Figures [5] a, b. 
Concluding remarks 

A brief overview over some recent results on the spectra of dynamo operators has been given. The obtained 
structural features like the resonance effects in the unfolding of diabolical points as well as the unexpected link 
to KdV soliton potentials, elliptic integrals, SUSYQM and the Dirac equation appear capable to open new 
semi-analytical approaches to the study of a 2 — dynamos. 
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